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Abstract 

We give a set of exact nonlinear closed-form solutions for the non-spherical col¬ 
lapse of pressure-less matter in Newtonian gravity, and indicate their possible 
cosmological applications. 
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1 Introduction 

The Zcldovich approximation [1], of decades ago, serves as a reliable basis for mod¬ 
ern studies of the large-scale cosmic structure formation (see the classic paper [2], 
recent works [3], and the references therein). It stemmed from an elegant exact 
nonlinear solution for the Newtonian anisotropic motion of dust-like gravitating 
matter, which describes a growth of initially weak adiabatic perturbations due to 
gravitational instability in an expanding universe. 

We re-visit the Zeldovich solution and its generalizations mn to reveal their 
new features emerging under the time-reversal transformation (TrT), t —y —t. To 
start with, we note that Friedmann’s general-relativistic cosmological solution is 
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known to have a Newtonian analog corresponding to a uniform isotropic expansion 
of a perfect pressure-less gravitating fluid (‘dust’). After the TrT, it turns to the 
solution describing the isotropic gravitational collapse (parabolic motion): 

x{t, 0 = {t 0 - f) 2/3 £ p(t) = [6tt G(t 0 - t) 2 ] -1 . (1) 

Here x is the Euler, and £ is the Lagrange coordinate, G is the Newtonian gravita¬ 
tion constant, p(t) is the matter density, and to is the collapse time, when the scale 
factor a(t ) = (t 0 — t) 2 / 3 becomes zero and the density goes to infinity. Solution ([I]) 
applies to uniform spheres of a finite total mass and any initial radius. 

In this paper we give a set of solutions for the non-spherical (anisotropic) col¬ 
lapse using the time reversal of expansion solutions found in [TJ SJ |5], and discuss 
their properties. The main reason for this is that in modern cosmology three types 
of building blocks are recognized forming the large-scale Cosmic Web. Those are 
rich spherical clusters of galaxies, flat superclusters (‘pan-cakes’), and elongated 
filaments. With all the necessary reservations, the formation and evolution of 
spherical clusters can hopefully be described by the asymptotically spherical col¬ 
lapse solution (HD • Superclusters and filaments are seen as highly oblate and highly 
prolate structures, respectively. Their evolution may be described by the solutions 
o.o and containing the growing mode, assuming that they correctly 

expose the basic features of a real finite mass collapse with the same dynamical 
asymptotics. Generalized formulas (fT5D add dark energy to the picture. 

One can construct the collapse solutions directly, without the TrT; this explains 
the powers of t = to — t involved in them. Assume that the flow (JT]) is perturbed 
in one direction by a term with separated variables: 

xi = (t) 2/3 £i + <7 (t)/(£i); £ 2i3 = (a) 2/3 &,3 ; (2) 

functions g(r) and /(£i) are so far arbitrary. By repeatedly differentiating these 
expressions one obtains the accelerations d 2 x n /dt 2 1 n = 1,2,3; according to the 
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equations of motion, they are the negative of the respective components of the grav¬ 
itational potential gradient. The divergence of the gradient produces the Laplacian 
of the potential; divided by AttG, it gives, by the Newton law of gravity, the first, 
‘gravitational’, formula for the matter density. The other one comes from the mass 
conservation, i.e., the continuity equation; the two densities coincide when 

9 ~ (4/9 r 2 )g = 0 ; (3) 

linearly independent solutions of this equation are r 4//3 and r^ 1 / 3 . The first one 
corresponds to the solution ((U) of sec. [21 the second gives the solution ([2]), sec. El 
A linear combination, t 4 / 3 /(£i) + t _1//3 F(£i), turns out also possible in the expres¬ 
sion (J2|) for xi, allowing for the solution (ITU]) of sec. [4] 

2 Collapse with decaying asphericity 

The well-known Zeldovich’s ‘pan-cake’ solution [1] for planar motion of pressure- 
free matter takes the form, after the TrT : 

xi(t, £i) = (*o - t) 2/3 £i + (to ~ C /3 /(£i), x 2,3(t, £2,3) = (to -1) 2/3 6,3 ; (4) 

^ = 6nG(to - t ) 2 l + (t 0 -t) 2 / 3 /'(6)’ f ^ = d£i' (5) 

The deviation in x\ from the spherically symmetric motion (JT|) goes to zero when 
t —$■ t 0 — 0, becoming negligible against the background (a decaying mode). So 
formulas (]T|) . (l5|) describe an initially non-spherical collapse which undergoes spher- 
ization and becomes completely spherical in the end. 

The density (JUJ) is non-singular until the collapse when 1 + (t 0 — t) 2 t 3 f(£i) > 0 
for all £ 1; guaranteeing also a unique inverse, £ = £(i, x), of the law of motion 
x = x(t, £), as required. The inequality, and thus the solution, is apparently valid 
for a finite period of time, U < t < to, if /(£ 1 ) > —(to — ti)~ 2 t 3 . For = —00, 
function /(£ 1 ) is non-decreasing, / (£i) > 0, so the solution holds on the whole 
semi-axis —oo<t<to- 
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3 Collapse with growing asphericity 

A counterpart to the solution (H|) . (l5]l comes from |[3]: using the TrT, one finds: 


xi (t, £1) = (to ~ t) 2/3 £i + (t 0 - t) 1 / 3 A(£i), x 2 , 3 (t, £2,3) = (to - t ) 2/3 £ 2 ,3 ; ( 6 ) 

p(i ’ 6) = feG0^^rT(i^pFK0 ■ (7) 

The deviation from the spherical symmetry here increases with the time (grow¬ 
ing mode); initially (t = — 00 ) the flow is entirely spherical. The density is 
non-singular until t = to if and only if F'(£ 1 ) > 0. Unlike the previous two 
cases o and near the collapse it is inversely proportional to the first, in¬ 
stead of the second, power of (to — t) and, generically, remains non-uniform, 
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6ttG( t 0 — t)F'(£ 1 ) ; instead of a point, the dust collapses to the oq-axis. 

If the last inequality is violated, and the minimum of F'(£ 1 ) is at a single point, 


mmF\^) = F\Ci) = -F<0, ( 8 ) 

a 

then a local collapse occurs at t = t 0 — t* before the global one at t = t Q : the 
density becomes infinite at the plane £1 = £*, or aq = t 2 J 3 ^l + U 1 ^ 3 F(£J : ). Note 
that p(t, ££) is maximum before the collapse; at the collapse, t = t 0 — t*, the density 
at all other points £1 7 ^ £ 5 " remains finite and non-zero. 

If the same minimum (jSJ) is achieved at a number of discrete £i-points, then 
the ‘early’ collapse happens at all the corresponding planes. If Eq. (jSJ) holds for 
some £ 1 - interval, 

F(£0 = —t*£r, a < £1 < P , (9) 

then xi(t, £ 1 ) = (t 0 - t) _ 1 / 3 (t 0 - F - t)£ 1 , x 1 (t 0 - t*,£ 1 ) = 0 for a < £1 < j3, 
and this whole £1 -layer collapses to a single plane aq = 0. Ultimately, when 
Eq. (jSJ) is valid for all £1, i.e., F(£i) = —t*£i, the early collapse becomes global. 
the density p(t, £ 1 ) = [67rG(to — t)(t 0 — t* — £)] -1 is uniform, so the whole space 
collapses. Still, the flow remains anisotropic and the dust goes to the plane aq = 0, 
and not to a point. 
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4 General solution 


In paper [3], a solution containing both modes was also found; its time reversal is: 
x\ (t, Ci) = (^o - t) 2/3 Ci + (t 0 - t) 4/3 /(Ci) + (t 0 - t) _1/3 F(Ci) ; (10) 

= 6iG(i„ - i) 2 1 + (i„ - Wm i) + («o - 0-V'Ki) (11) 

(x 2 .3 are as before). It is general in a sense that the two free functions /(£i) and 
F(Ci) allow one to meet any initial conditions for both the Xi-coordinate and X\- 
velocity. The solution (TTOji . (TTT]i is valid at all times for f'(£i) > 0, F'(£ i) > 0; 
other cases are analyzed as above. Of course, the decaying mode plays no role 
near the collapse, which is thus either an early (t < f 0 ) planar, or a global axial 
one at t — to- 

Note that solutions describing the expansion with an anisotropic flow and uni¬ 
form density p(t) = [67 t G(t — to)] 2 were studied in paper [6]. Their time reversals 
can add to the set of non-spherical collapse solutions obtained here. 


5 Anisotropic collapse of a finite mass 

All the solutions of secs. [2] - [4] apply to circular cylinders of a finite mass with the 
symmetry axis along x\, of an arbitrary initial radius and height, whose density is 
uniform in both the radial and azimuthal directions. Such a cylinder is specified 
by the following bounded set of Lagrange variables: 

a < Ci < /3, i/cl + Cl < 7; -°o < a < (3 < oo, 0<7 <oo. 

The decaying-mode solution of sec. [^describes the collapse of a cylinder to a point. 
Two other solutions, of secs. [3] and [4J, containing the growing mode, end up with an 
early (t < to) collapse to a disk (or a number of disks) perpendicular to the x± axis, 
when F'(£i) < 0 for some values of ^ in the above range. Otherwise ^ (^i) > 0 
for all £i, and the cylinder collapses at t — to to a part of the xi-axis. Depending 
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on the behavior of F{£i), it can be either the whole axis, or semi-axis, or a finite 
segment moving to infinity when t —> f 0 — 0. So the outcome of the collapse process 
is determined by the initial conditions in x \ -direction, if other initial conditions 
provide for the axial symmetry. Namely, at an initial moment of time t — ti the 
velocities should be consistent with the coordinates, v 2;3 = —(2/3)[2:2,3/ (t 0 — i*)], 
which relations hold then on all the way to the collapse. 


6 Collapse on the dark energy background 


The original Friedmann cosmological solution includes the cosmological constant 
A, which may be zero or non-zero. The time-reversed isotropic flow with A > 0 is: 


x(t,£) = a(t)i, pit) oc [ait)} 3 ; 


a(t) oc sinh 2//3 |^(3/2)(\/3A)(t 0 — t) 


( 12 ) 


The corresponding non-spherical solution of paper [5] after the TrT becomes: 


xiit, 6) = a(*)fi + 1) + &(t)F{£ 1), ar 2 , 3 (t, 6,3) = a(t)6,3 • ( 13 ) 

Functions 0(t) and $(t) giving the time dependence of the decaying and growing 
mode are found in [SJ. Solution (1T5D is general in the same sense as in the case (TT01) . 
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